Splitting of doubly quantized vortices in dilute Bose-Einstein condensates 
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We investigate the dynamics of doubly charged vortices generated in dilute Bose-Einstein conden- 
sates by using the topological phase imprinting technique. We find splitting times of such vortices 
and show that thermal atoms are responsible for their decay. 



In their recent Letter Huhtamaki et al. theoret- 
ically investigated the splitting of a topologically im- 
printed doubly charged vortex into two singly charged 
vortices as occurring in a dilute atomic Bose-Einstein 
condensate. They compare the results of simulation with 
recent experiment [|J and show that the combination of 
gravitational sag and the time dependence of the trap- 
ping potential alone are enough to explain the observed 
splitting times. Based on such an outcome the authors 
of Ref. [l| claim that, contrary to previous theoretical 
results [3J, the thermal excitations are not relevant in 
modeling the experiment of Ref. Q. We are going to 
show in this Brief Report that, indeed, the opposite is 
true. In fact, a number of thermal (uncondensed) atoms 
appears in the system while disturbing the gas. They 
continue to appear after the perturbation is over and un- 
til the doubly quantized vortex breaks into two singly 
quantized vortices. The overall number of uncondensed 
atoms remains approximately on the level of 20%, which 
is already at the edge of experimental detection capabil- 
ities. However, the uncondensed atoms do not form the 
broad cloud allowing the identification by fitting to a bi- 
modal distribution - they are rather located in the core of 
the vortex and therefore are harder to detect. Perhaps, 
the signatures of the presence of thermal atoms in vor- 
tices cores are already visible in experiment in a way that 
after splitting the cores of two singly charged vortices get 
darker in comparison with the core of initially imprinted 
doubly quantized vortex (see Fig. 2 in Ref. 0]). 

To investigate the thermal excitations in a Bose gas we 
use the classical fields approximation [J - an approach 
that treats both condensed and thermal atoms at the 
same footing until the detection time when the split- 
ting into the condensate and the thermal cloud occurs. 
Technically, such a decomposition requires calculation of 
time and space average of a one-particle density matrix 
built of the classical field evolving according to the Gross- 
Pitaevskii equation [1]. 

Therefore, we have repeated the calculations of Ref. 
Q. As in (lj, we closely follow the experiment reported 
in 0] and solve the Gross-Pitaevskii equation with time- 
dependent trapping potential (according to the idea of 
topological phase imprinting ||) combined with the grav- 
itational potential. However, we interpret the solution of 
the Gross-Pitaevskii equation as a classical field and by 
using the space averaging procedure we determine the 
condensate and the thermal cloud at each time. 



FIG. 1: (Color online). The z-integrated (over [—15 /im, 15 
fj,m] interval) condensate (a), thermal (b), and total (c) den- 
sities at 47 ms for the interaction strength an z = 8.5; the 
horizontal cuts of the condensate (dashed line), thermal (dot- 
ted line), and total (solid line) densities (d); the condensate 
fraction as a function of time (e). 



Fig. Q] clearly shows that the uncondensed atoms ap- 
pear in the system during the evolution. Although ini- 
tially all atoms are in the condensate (i.e., the condensate 
fraction equals 1 as in Fig. [He)), already after 6ms dis- 
tinguishable fraction of uncondensed atoms is produced. 
This is not surprising since the process of imprinting the 
vortex is accompanied by a sudden squeeze of the con- 
densate in the radial direction and a kick of it in the 
vertical direction Q . The thermal atoms continue to ap- 
pear after the imprinting is over and until the doubly 
quantized vortex splits into two singly charged vortices 
(production of thermal atoms while the system was ini- 
tially at zero temperature was also reported in Ref. [||, 
where the crystallization of a vortex lattice in a stirred 
Bose-Einstein condensate was investigated). 

Fig. [U also proves that the thermal atoms are mainly 
located in the vortex core. On the other hand, the ther- 
mal noise on the level of about 20% leads to the decay 
of a vortex in approximately 45 ms (see Fig. 3 in Ref. 
Q) in agreement with the calculations in [l[ and this 
Brief Report. We claim this means that, indeed, be- 
hind the origin of the splitting of the doubly quantized 
vortex reported in [l[ (the impetus given by the time de- 
pendence of the external potential) lies the presence of 
uncondensed atoms appearing in the system as a result 
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of its disturbance. In other words, the understanding 
of the experiment by Shin et al. requires going beyond 
the standard Gross-Pitaevskii approximation. It is clear 
that other kinds of perturbation of the condensate (for 
instance as the one considered in Ref. |7]) will result in a 
production of uncondensed atoms and consequently lead 
to the decay of the multiply charged vortex. 

In conclusion, we have studied the splitting process of 
doubly charged vortices created via the topological phase 
imprinting method. We found that the uncondensed 



atoms, inevitably produced in the condensate while it 
is disturbed, are responsible for the decay of such vor- 
tices. This is because the unstable modes, located in the 
vortex core, are seeded by the uncondensed atoms which 
are also located there. 
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